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1. Introduction 



Light-by-light scattering is one of the most fundamental processes in QED. Theoreti- 
cally, it proceeds at leading order, C(a 4 ), via one-loop box diagrams containing charged 
particles. At center-of-mass energies \fs far below the mass of the electron, the process 
is described by the Euler-Heisenberg effective Lagrangian and the cross section rises 
rapidly with energy, a oc s 3 /m*. The cross section peaks at ~ 3m e 0, then begins 
to fall rapidly, a oc 1/s at fixed angles for s ^> ml || . At still higher energies, similar 
thresholds are crossed for the muon, tau, and light quarks - or rather, the light hadrons. 
(See fig. |l|.) The final significant Standard Model thresholds reached are those of the 
W boson (|, H and the top quark. 

The direct experimental evidence for 77 — > 77 scattering is still scant, particularly 
for energies above a few GeV. At optical (electron volt) energies, at least two experi- 
ments have been proposed to detect the Euler-Heisenberg interaction via birefringence 
of the vacuum using laser photons in a magnetic field || . The Crystal Ball experiment 
at the SPEAR storage ring had some unpublished evidence for 77 — > 77 scattering 
at y^i of order several MeV, after subtracting beam-off 7 ray backgrounds 0. More 
extensive data are available for the process of Delbriick scattering, in which two of the 
four photons are supplied by the Coulomb field of a nucleus. Several experiments have 
studied this process, for incident photon energies ranging from a few MeV to about 
a GeV. Higher order Coulomb (Za) corrections are often important in the compari- 
son with theory §|. Light-by-light scattering via an electron loop is also tested quite 
precisely, if indirectly, by the measurement of the anomalous magnetic moment of the 
electron |§, as well as that of the muon PI (the latter is even sensitive to muonic and 
hadronic light-by-light scattering). 

At much higher energies, the small size of the Standard Model light-by-light scat- 
tering cross section provides a potential window to new physics. By backscattering a 
laser pulse off an intense, high energy electron beam fLl| , it is possible to create 77 
collisions with yfs of order 100-1000 GeV, high levels of initial state polarization, and 
luminosities of order tens of fb _1 per year. It has been shown that the 77 — > 77 process 
at high energies is a sensitive probe |12|] of theories with large extra dimensions |13 



for example. Despite this interest in light-by-light scattering, the process has to date 
only been calculated to leading order (one loop), except in the low energy limit, where 
the two-loop corrections to the Euler-Heisenberg Lagrangian are known JT4 ] . 

In this paper, we present the two-loop corrections to 77 — > 77 in the ultra- 
relativistic regime where the kinematic invariants s, t, u are much greater than the 
(squared) charged fermion masses. This regime is relevant for two ranges of center of 
mass energy in the Standard Model: 
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Figure 1: The leading order (one loop) cross section for light-by-light scattering, for inter- 
mediate energies, showing the light fermion thresholds. The cross section has been integrated 
over a range of center-of-mass scattering angles, 30° < 9 < 150°, for the two independent 
choices of initial photon helicities (assuming CP invariance), ++ (solid) and H — (dashed). 
We have multiplied a by s merely to compress the dynamic range. We set a = 1/137.036. 
For yfs < 0.5 GeV we omit all quark loops but include the charged pion and kaon loops. For 
y/s > 0.5 GeV we use the quark loops, with m u = = 5 MeV, m s = 100 MeV, m c = 1.25 
GeV, m b = 4.2 GeV, m t = 175 GeV, and m w = 80.42 GeV. The hadronic contribution to 
the region from 0.5 to 8 GeV is not computed reliably by the quark boxes; hence that region 
is shown dotted. We also omit the hadronic resonance contributions, 7r°, rj, etc. 



• For m c <C < 2m^ (and neglecting the tiny bottom quark contribution), the 
QCD corrections from attaching a gluon line to the quark box are most important. 
These give rise to finite 0(a 4 a s ) corrections to the cross section. 

• For m e <^ y/s < 2m ^ 2m n , the dominant corrections are QED corrections from 
attaching a photon line to the electron box. There are also QED corrections 
from inserting a fermion loop onto an external leg; however, these are cancelled 
completely if the theory is renormalized at zero momentum transfer, i.e. by using 
a = a(0) = 1/137.036... as the coupling constant. Then the two-loop QED 
corrections become identical to those of QCD, up to an overall constant. 
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The Feynman diagrams for 77 — > 77 at two loops are a small subset of those 
required for gluon-gluon scattering, gg — > gg. The interference of these two loop am- 
plitudes with the gg — > gg tree amplitudes is an essential ingredient for obtaining the 
next-to-next-to-leading order QCD corrections to jet production at hadron colliders. 
This interference was recently evaluated in a tour de force calculation by Glover, Oleari 
and Tejeda-Yeomans [r5| . In the light-by- light case, the tree amplitudes vanish. Thus 
the nexi-to-leading order corrections require a different interference, of two-loop am- 
plitudes with one-loop amplitudes. Instead of evaluating this interference directly, we 
have computed the two-loop 77 — > 77 amplitudes in a helicity basis. Thus polarized as 
well as unpolarized cross sections can be obtained, information which is quite useful in 
the photon case, because of the high initial state photon polarizations that are possible 
with backscattering. 

Rather than calculating the Feynman diagrams for the two-loop 77 — > 77 helicity 



amplitudes, we employed a unitarity- or cut-based technique [T^, [T7|, to generate 
the required loop momentum integrals. These integrals were then evaluated using 



techniques recently developed |2(J, £1]] to handle double box and related integrals 
where all internal lines are massless. The restriction to massless internal lines limits 
the validity of our results to the ultrarelativistic region where all kinematic invariants 
are much greater than the relevant fermion masses. 

Remarkably, some of the helicity amplitudes stay quite simple, even at two loops. 
Some of this simplicity can be understood via unitarity and supersymmetry Ward 
identities |22| . 

This paper is organized as follows. In section 2 we present analytic results for 
the two-loop 77 — > 77 helicity amplitudes. (The more complicated formulae are rele- 
gated to an appendix.) In section 3 we give numerical results for the QCD and QED 
corrections. In section 4 we present our conclusions. 



2. The Two-Loop Amplitudes 

The one-loop Feynman diagrams for photon-photon scattering via a charged fermion 
consist of six box diagrams, which are related by permutations of the external photons 
to the diagram shown in fig. |2|a. At two loops, there are 6 • 10 = 60 nonvanishing 
Feynman diagrams, corresponding to connecting any two sides of the box with a photon 
propagator (or a gluon propagator, in the case that the charged fermion is a quark). 
Up to permutations, these diagrams are depicted in fig. ||b. Diagrams consisting of two 
separate charged fermion triangles connected by a single photon or gluon, as illustrated 
in fig. 0c, vanish by Furry's theorem |23] or simple group theory. 
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Figure 2: Feynman diagrams for light-by- light scattering via a charged fermion loop, (a) 
One of the six box diagrams contributing at one loop; the remaining diagrams are obtained by 
permuting the external photons, (b) The three types of diagrams contributing at two loops; 
the rest of the 60 diagrams are obtained by permutations. The curly line can be either a 
photon or a gluon, depending on whether the QED or QCD correction is being computed, (c) 
These diagrams vanish by Furry's theorem or simple group theory, (d) These diagrams could 
contribute to the QED correction. However, they are precisely cancelled by conventional 
on-shell renormalization. 



We did not evaluate the Feynman diagrams directly. Instead we computed the 
unitarity cuts in various channels, working to all orders in the dimensional regulariza- 
tion parameter e = (4 — D)/2. Essentially we followed the approach first employed at 
two loops for the special cases of iV = 4 supersymmetric amplitudes jl7|] and the pure 
gluon four-point amplitude with all plus helicities These amplitudes were simple 
enough that a compact expression for the integrand could be given. The fermion loop 
contributions with all plus helicities are about as simple ||24|| . However, for the generic 
helicity configuration, the integrands become quite complicated. 

We have used general integral reduction algorithms developed for the all-massless 
planar four-point topologies |19|, |20f , in order to reduce the loop integrals to a minimal 
basis of master integrals. Some mild extensions of these techniques are required in 
order to incorporate polarization vectors for photons of definite helicity f2f|. We then 
expand the master integrals in a Laurent series in e, which begins at order 1 / e 4 . Many of 
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the master integral Laurent expansions quoted in refs. []19, are in terms of Nielsen 



functions [p6fl , usually denoted by S n>p (x), with n + p < 4. However, using various 



identities |27| the results can be expressed in terms of the polylogarithms 

00 x i r x dt 



Li n (x) =E^= / -r Li n-i(i), 
i=1 1 jo 1 

Li 2 (*) = - r^ln(l-t), (2.1^ 



with n = 2,3,4. 

It is reassuring that all of the poles in e cancel for each helicity amplitude in the 
QCD case. In QCD, two loops is the first order at which a s appears in the 77 — > 77 
amplitude; therefore there can be no ultraviolet divergence. Any infrared divergence 
would have to be cancelled by real gluon radiation; but the process 77 — > 77$ is 
forbidden by group theory. 

In addition to infrared finiteness, the calculational framework was also tested by 
gauge invariance: replacing a photon polarization vector by its momentum vector pro- 
duces a vanishing result. Finally, the same computer programs for evaluating the cuts 
and reducing the integrals were also used to compute the fermion-loop contribution to 
the two-loop gg — > gg amplitudes in the helicity formalism and the 't Hooft-Veltman 
dimensional regularization scheme p^j . The interference of the two-loop gg — > gg he- 



licity amplitudes with the tree amplitudes, after summing over all external helicities 
and colors and accounting for the different scheme used, is in complete agreement with 
the calculation using conventional dimensional regularization ||15|| . 

The QED case requires exactly the same set of two-loop diagrams as in QCD, up to 
an overall factor. In addition, there are external fermion bubble insertions of the form 
shown in fig. §i. In dimensional regularization with massless fermions, these diagrams 
would vanish by virtue of containing scale-free integrals. This vanishing represents 
a cancellation of ultraviolet and infrared divergences. However, if one renormalizes 
QED in the conventional on-shell scheme, to avoid infrared divergences one should 
retain a fermion mass in the external bubbles. Now the bubble integral is nonzero and 
ultraviolet-divergent. But this divergence, and indeed the entire integral, is exactly 
cancelled by the on-shell-scheme counterterm, precisely because the external leg is a 
real, on-shell photon. 

Hence the on-shell-renormalized two-loop QED-corrected amplitude is the same 
as the two-loop QCD-corrected amplitude, up to overall coupling constant factors. In 
the on-shell scheme, the coupling constant should of course be set to a = a(0) — 
1/137.036.... This value should be used for all the QED couplings associated with 
the real, external photons — that is, for all couplings except the extra one associated 
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with the virtual photon in each two-loop graph. The typical virtuality of the extra 
photon is not zero, but of order s, assuming that the kinematic invariants s, t,u are 
all comparable in magnitude and much larger than the squared fermion mass, raj. 
Thus a running coupling a(/x) with \x ~ \fs should be used for the virtual photon 
insertion; whether it should be the MS running coupling, or that defined via the photon 
propagator at momentum transfer /i, is theoretically indistinguishable at this order. For 
m e ji < m M , the latter coupling at one-loop order is 



= r^rr-\ — rr- ( 2 - 2 ) 

[ln(^) " 



l - — 

3vr 



For n > Im^, the running of a(fi) receives hadronic corrections as well, which are best 
evaluated via a dispersion relation using the e + e~ — > hadrons data. In any event, 
the precise a(/x) used makes very little difference, since it only appears in an 0(a) 
correction to a process with a rather tiny cross section. 

For the QCD corrections, we use the MS running coupling a s (fj), again with ft = 
y/s, keeping in mind that our calculation is only to leading order in a s . 

We consider the process 

7(*i, Ai) +7(*2,A 2 ) -> 7(^3, A 3 )+ 7(^4, A 4 ), (2.3) 

where ki and A« are the photon momenta and helicities. In terms of the center-of- 
mass energy ^/s and scattering angle 9, the Mandelstam variables are s = (k\ + /C2) 2 , 
t = (k\ — k±) 2 = —s/2x (1-COS6 1 ), andw = (/ci-/c 3 ) 2 = -s/2x (1 + COS6 1 ), with s > 0, 
t < 0, u < 0. Parity, time-reversal invariance, and Bose symmetry imply that there 

are only four independent helicity amplitudes, M h+ , M_ +++ , M ++++ , and M H (.. 

Actually, crossing symmetry relates M ++++ and M_| |_. However, representing the 

two-loop amplitudes in a crossing-symmetric fashion in terms of master integrals would 
lead to more cumbersome formulae, so we shall present all four amplitudes instead. We 
adopt the overall phase convention of refs. 0, |29|. Then formula (9) in ref. |29| can 



be applied to our helicity amplitudes in order to obtain the 77 — > 77 differential cross 
section for generic circular and transverse photon polarizations (Stokes parameters). 

The one- loop helicity amplitudes due to a fermion of mass mj in the loop take a 
very simple form in the ultra-relativistic limit, {s,t,u} 3> m'j |3], |2], |3(|. We write 



M ^ = 8 iVQ 4 a 2 M« 2A3A4 , (2.4) 

where N is the fermion color factor (3 for quarks, 1 for leptons), and Q is the fermion 
charge in units of e. The functions are given by 

M« ++ = l, 
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M. 



(i) 



m2> lt2 + " 2 



2 s 



2 



ln^-Wl-^lnf-)-! 



t — U , / t 



Mil, = ^ln 2 ( — ] ln( — -1 



t 2 + s 2 / t \ t-s 



VK 



It 2 v s' u 



(2.5) 



The QCD- and QED-corrected two- loop amplitudes are 

M^Zir = 4 (iV 2 - 1)Q 4 a 2 ^ M{f A2A3A4 , 

*C3£r = 8 N Q 6 « 2 ' (2 - 6) 

where the explicit values of the M^ } x are given in Appendix ^|. 

Note from eq. ( |A.1| ) that M |_ + remains remarkably simple even at two loops — it 

is just a constant, independent of s, t and u\ This simplicity is actually predictable. The 
lack of an imaginary part for M^lL ++ (as for M^l1 ++ and Mi 1 | ++ ) can be deduced by 
considering the possible unitarity cuts of the amplitude. Let / be a massless fermion. 
Then a supersymmetry Ward Identity |22| shows that the tree-level amplitudes for 



77// and 777// vanish if all the photons have the same helicity (when considered as 

outgoing particles). The same result is true if one (or more) photons are replaced by 

gluons with the same helicity. Assigning opposite (outgoing) helicities to intermediate 

particles on opposite sides of a cut, it is easy to see that all two- or three-particle cuts 

(2) 

vanish for the all-outgoing-plus-helicity two-loop amplitude M__ ++ . This argument 
assumes that a D = 4 helicity assignment is valid for the intermediate particles, a 
condition which is justified for the processes considered here by the absence of infrared 
or ultraviolet divergences. Thus M__ ++ must be a dimensionless rational function. 

Crossing symmetry implies that it is totally symmetric in s, t, u. Imposing at most 

(2) 

single poles in s, t, u, and using s + t + u = 0, one finds that M__ ++ must actually be 
a constant. 

Equation (|A.2 ) shows that Mi 2 + ++ is also fairly simple, containing only logarithms, 



and not Li2, Li3 or Li4 functions. Its two-particle cuts in D = 4 are the product of 
two rational functions (just like the cuts of a one- loop amplitude), which may partially 
account for this simplicity; however, the three-particle cuts do not similarly simplify at 
a glance. 

Interestingly, up to the overall normalization these amplitudes are identical to 
the subleading color contributions for the gg — > 77 amplitudes presented in ref. fl31| . 
(In that paper all particles are taken to be outgoing, which corresponds to flipping 
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the helicity labels for particles 1 and 2.) These amplitudes are relevant for improved 
estimates of the di-photon background to production of a light (mass < 140 GeV) Higgs 
boson at the Large Hadron Collider ||32|| . 

It is instructive to quote the values of the helicity amplitudes in various limits. The 
values of the one-loop functions for 90° scattering (t — —s/2, u = —s/2) are 

M a 2 ++ (90°) = 1, 
M« ++ (90°) = 1, 
Mj+ ++ (90°) « -3.46740, 

M| 1 2_ + (90 ) « -0.12169 + 0.46574vri , (2.7) 

while the two-loop functions are 

Mi 2 ! ++ (90°) = -1.5, 

M {2) +++ (90°) « 0.17770 - 0.23287m , 

MS ++ (90 o ) w 1.24077+ 1.00717tu, 

M| 2 2_ + (90 ) « 0.01445 + 0.36840vri . (2.8) 

Note that except for the first, the two-loop functions are smaller than their one-loop 
counterparts. 

In the small-angle limit \t\ <C s, the values of the one-loop amplitudes, up to 
corrections suppressed by powers of t/s, are 

m! 1 ! ++ ~ i , 

Mi 1 ) ~ 1 , 
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~ -lx* - X - — -1, 
++++ 2 2 

Af£L+ ~ - l -X 2 - X - 1 - m(X + 1) , 

a4 1 2 + _ ~ , 

(2.9) 

while the two-loop amplitudes are, 



r(2) _ , 3 



2 ' 

M (2 | ++ ~ -X' + -X + - - - + m ( -X + - , . 
+++4 2 8 4 V4 4 1 

r 2 



24 12 Vi2 2> V12 2' 360 s 12 4 
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+ m (Ix 2 + ix + 2C 3 + — - 



M. 



(2) 



M 



(2) 



r\j 



lx 4 - -X 3 
24 12 

+ VK 

1 

2 



(—X 3 



(- 

V 6 



4 



X 2 - 

*X 
2 



13 

C) T 

5 „ 



-7T 



6 ' 2 
2C 3 + ^ 



X 



11 
180' 



-7T 



5 2 1 

-7T 



12 



(2.10) 



where X = — ln(— s/t). Notice that the small-angle scattering amplitudes where both 
photons flip their helicity have no logarithmic enhancements at one- or two-loops; those 
with one helicity flip have no logs at one loop but a factor of X 2 at two loops; and those 
with no helicity flips pick up one factor of X 2 for each loop. For small angles these 
contributions are actually power suppressed compared to t-channel vector exchange 
diagrams which first appear at three loops p3| . 

Next we give the differential cross sections in terms of the above amplitudes. The 
leading-order 77 — > 77 unpolarized differential cross section for a single fermion flavor 
is given by 



da LO 
dcos9 

where 



N 2 Q* 



a 



2ns 



\M. 



(i) 



+ 4IM. 



(i) 



+ \M. 



(i) 



+ \M_ 



(i) 



\M. 



(i) 

+-+- 



(2.11) 



M { + L 2 + _{s, t, u) = M| L) _ + (s, u, t) . (2.12) 

The QCD- and QED-corrected unpolarized differential cross sections for a single ferm- 
ion flavor are given by 

da QCD da LO 



dcos9 



+ 



da a ° 



d cos 9 d cos 9 
da LO da a 

d cos 9 d cos 9 d cos 9 



(2.13) 



where 



da c 



dcos9 
da a 

dcos9 



x(x 2 -i)g 8 ^^/( 2 - 

27H S 

N 2QlQ gVri j(a,l) 
7T 2 S 



1) 



(2.14) 



and 



1^ = Re 



M. 



(2) 



M. 



(i) * 



4M. 



(2) 



M. 



(i) * 



M. 



(2) 



M {1) 



+ M. 



(2) 



M. 



(i) * 



M. 



(2) 



M. 



(i) * 



(2.15) 
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The cross sections for circularly polarized photons can be constructed easily from 
eq. ( p. 14j ) by reweighting the products of M's. For arbitrary initial photon polar- 
izations, linear as well as circular, eqs. (9)-(15) of ref. can be used, after making 
the replacement (for the QCD case) 

F Xl x 2 x 3 x 4 - -M Ai a 2 a 3 a 4 = 8 NQ 4 a 2 + 4 (N 2 - 1)Q 4 a 2 ^ M$?^ ■ 

(2.16) 

(Note that our definitions of t and u are reversed with respect to ref. [29].) 



3. Numerical Results 

The QCD K factor is conventionally defined as the ratio of next-to-leading to leading- 
order cross sections, 

da*™ /d cos 9 
~ da LO /d cos 9 ' 1 ' 

In a region between quark thresholds, and below the W mass, we can write 

g = 1 + s^ aw-i)q? fc , (3 . 2) 

Ei NiQi 



where 

E{A} Re 



k 



M (2) M (1) 



E W |M { (1) ' 2 



(3.3) 



<{A}I 

The label % runs over the number of quarks and leptons with masses much less than 
yfs\ N{ = 3 for quarks, Ni = 1 for leptons; and Qi is the fermion charge. For energies 
between m& and mw, for example, (Yli(N? — l)Qf) / (Ysi NiQf) = 70/87. For energies 
above mw, one should add the W loop contributions to Ai 1 ^ ^. This dilutes the QCD 
corrections considerably, since the W loops quickly dominate 0, ||] (see fig. [l|), and 
they have no QCD corrections. 

In fig. |3|, k is plotted as a function of the center-of-mass scattering angle 9. Due to 
Bose symmetry, only the forward region cos 9 > has to be plotted. We also plot the 
corresponding curves for fully circularly polarized initial photons, but summing over 
final-state helicities, 

Ea 3 ,a 4 R- e 



&AiA 2 — 



M (2) M (1) 
m { X y vl { X } 



Ea 3 ,a 4 |M« 



2 



(3.4) 



Taking into account CP invariance, there are two independent cases, ++ and H — . 
Although the QCD corrections are of order a s /ir for these two cases, they have opposite 
signs, and in the unpolarized cross section there is a large cancellation. The QCD 



10 




Figure 3: Eqs. (]3.2|)-( |iO| ) express the QCD-correction K factor in terms of the quantities k 

(unpolarized initial and final photons), k ++ and fc_| (fully polarized initial photons). Here 

k (solid), k ++ (dotted) and (dashed) are plotted as a function of cos#. 



corrections to the unpolarized cross section are only of order 0.1 x a s /ir for central 
scattering angles, becoming as large as a s /ir only for | cos^l > 0.95. 

Figure |] displays the QCD corrections to the cross section in the 10-100 GeV region 
where quark boxes are important, and the quarks can be taken to be approximately 
massless. The QED corrections can be neglected in this region as they are a factor 
of 8 or so smaller. The two independent initial helicity configurations are shown, 
after integrating the differential cross section with an acceptance cut of 30° < 9 < 
150°, and multiplying by s as in fig. [l|. The QCD corrections were computed by 
multiplying the leading order one- loop result, which contains the full fermion (and 
W) mass dependence, by the K factor formula ( |3.2| ) evaluated for massless fermions 
in both numerator and denominator (which are of course integrated separately over 
6). We can apply eq. ( |3.2| ) even across the b quark threshold, simply because the b 
quark's small charge of —1/3 lends it a tiny contribution to the 77 — > 77 cross section. 
(Nevertheless, we interpolate (EiW - Wt) / &i N iQt) from 272/345 » 0.788 to 
70/87 ~ 0.805 across the threshold.) We see that the QCD corrections are numerically 
rather small for central scattering angles in this entire energy range. The small size of 
the correction explicitly demonstrates the reliability of the leading order prediction. 
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QCD Corrected 77 -> 77 Cross Section 




s 1/2 (GeV) 



Figure 4: QCD corrections to 77 — > 77 in the intermediate energy region 8 GeV < y/s < 
100 GeV where quark boxes are important, for a s (mz) = 0.118, and renormalization scale 
fi = -y/s. Plotted are the results at leading order (dashed), and including the QCD corrections 
(solid), for the same quark masses used in fig. [j]. 



The QED corrections can also be extracted from fig. |3|. In the region between m e 
and m^, for example, the corrections to the unpolarized cross section are 

da^ ED _ da LO 
d cos 9 d cos 9 

For the same angular cuts (30° < 9 < 150°) as in the QCD case, this amounts to only 
a 0.07% decrease in the cross section for ++ initial helicities, and a 0.35% increase for 
H — initial helicities. It will be a true challenge to measure the 77 — > 77 reaction to 
this level of precision. 



1 + 2 



71 



(3.5) 



4. Conclusions 

In this paper we presented the two-loop QCD and QED corrections to light-by-light 
scattering by fermion loops in the ultrarelativistic limit where all kinematic invariants 
are much greater than the relevant fermion masses. These corrections reliably give the 
leading Standard Model corrections for most energies below 100 GeV. The corrections 
are quite small numerically, showing that the leading order computations are robust. 
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Some of the helicity amplitudes remain quite simple even at two loops, as a conse- 
quence of unitarity and a supersymmetry Ward identity for tree amplitudes. 

To extend our results to regions where the kinematic invariants are comparable 
to the masses in the loops, the technology for computing two-loop double box inte- 
grals should first be extended to include massive internal lines, which seems feasible. 
Probably the most important application would then be to compute the electroweak 
corrections to the W box contribution to 77 — > 77, since that contribution dominates 
at high energies, where new physics contributions are most likely to be found. 
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A. Two-loop helicity Amplitudes 



The explicit expressions for the two-loop amplitudes appearing in eq. ( |2.6|) are 

r(2) = _3 
' " 2' 



Ml 



M. 



(2) 



1 



^4±{{X + mf + vr 2 ) + \{x 2 + y 2 )((X - Yf + vr 2 ) 

y 2 



--x)(X 



VK 



jt <-y u\ 



(A.I) 



(A.2) 



M. 



(2) 



-2x 2 



Li 4 (-x) + Li 4 (-y) - (X + m)(Li 3 (-x) + Li 3 (-y)) 



1 



7T 



+— X 4 - -X A Y + — XY vr 4 

12 3 12 90 

[x - y)[Lu(-x/y) - — Li 2 (-x 



+ i-X(x 2 -3XY + 



— X 



2Li 3 (-x) - Li 3 (-x/y) - 3C 3 - 2(X + m)Li 2 (-x) 
+ (X - y)(Li 2 (-x/y) + X 2 ) + ±-(b(X -Y) + 18i7r)((X - F) 2 + tt 2 



12 



-X(X 2 + 7T 2 ) -Z7T(r 2 + 7T 2 



+ - 



1 1 - 2a; 2 



7T 



4 y 2 



((X + ztt) 2 + vr 2 ) - ~(2xy + 3)((X - Y) + tt 2 ) + — 



12 



(A.3) 
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M 



(2) 



U 4 (-x/y) - Li 4 (-y) + -(X -2Y- m)(Li 3 (-x) - Q 



— (X 4 + 2M 3 - AXY 3 + F 4 + 2tt 2 F 2 ) + — 7T- 
24 v 7 360 



Li 4 (-x) - (4 - -(X + m)(Li 3 (-x) - Cs) 



7T 



+ — (Li 2 (-x) 



7T 



^X< 



48' 



-X 4 



x 



+ (2--1 

1/ 

f 2(2- + l 

y 



U 3 (-x) - (X + m)Li 2 (-x) + Cs - - \{X + F) 
Li 3 (-y) + (Y + z7r)Li 2 (-x) - Cs + ^(2F 2 + vr 2 ) 



--X 2 (X + 3m) 

o 



-i(2x 2 -, 2 )((X-r) 2 + 7 r 2 



vr 2 )- 



2-jf 

4a; 2 



7T 



((r + m)^ + vr 2 ) + — 



6 



+ -(2x + y 2 ) 



-(X + Z7r) + — (y + 27r) 

2/ « 



Here 



(A.4) 
(A.5) 
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